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f*^ ■ Gravitational collapse in (n+2) dimensional quasi-spherical space-time is studied for a fluid with 

f***) ' non vanishing radial pressure. An exact analytic solution is obtained (ignoring the arbitrary in- 

\ tegration function) for the equation of state p r = (7 — l)p. The singularity is studied locally by 

comparing the time of formation of apparent horizon and the central shell focusing singularity while 
the global nature of the final fate of collapse is characterized by the existence of radial null geodesic. 
It is revealed that the end state of collapse for D dimension with equation of state p = — p, for 
(D— 1) dimensional dust and (D— 2)dimension with equation of state p = p are identical. 
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I. INTRODUCTION 



5-H , 

Gravitational collapse is an important and challenging issue in Einstein gravity, particularly after the 
formation of famous singularity theorems [1] and cosmic censorship conjecture (CCC)[2]. Also from the 
perspective of black hole physics and its astrophysical implications, it is interesting to know the final outcome 
of gravitational collapse [3] in the background of general relativity. The singularity theorems only tells us about 

t-H ' the generic property of space-times in classical general relativity, but it cannot provide us about the detailed 
features of the singularities whether an external observer can visualize the singularity or not. Moreover, the 
CCC is incomplete [4,5] in the sense that there is no formal proof of it in one hand and on the otherhand there 
are counter examples of it. It should be noted that choice of initial data [6] has a role in characterizing the 

Q\ • final state of collapse. 

o: 

Though a lot of works have been done for dust collapse [7-12] both for spherical and quasi-spherical models [13- 
18] in four and higher dimensional space-times but there is not much progress in studying gravitational collapse 
^ . in quasi-spherical space for perfect fluid or matter with anisotropic pressure [19]. In recent past Dadhich etal 
[20] examined the role of the equation of state in charecterizing the final state of spherical collapse and showed 
a similarity among collapsing processes in different dimensions with different eq. of state. In this work we have 
■ studied quasi-spherically symmetric collapse of a fluid with non-vanishing radial pressure (with eq. of state 
p r = (7 — l)p) in higher dimensional space-time. 

II. HIGHER DIMENSIONAL SZEKERES' MODEL 

The quasi-spherical Szekeres' space-time in (n+2)dimension has the metric ansatz (in the comoving coordi- 
nates) 

ds 2 = dt 2 -e 2a dr 2 -e 2 ^dx 2 (1) 

1=1 

where a and (3 being functions of all space-time co-ordinates are of the form 

e a = R' + Ru', e /3 = Re l/ with R = R(r,t) and v = v(r,xi, ...,x n ) (2) 
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For the matterfield with anisotropic pressure i.e., T£ = diag(p, — p r , — p T , 
equations (after some manipulation) 



., —p T ), we have from the Einstein 



Pr = - 1 h 



11 = &■ 



(3) 



F{r,t) = ^-R n - l e^ +1 >{R 2 - f{r)) 



(4) 



while from conservation equation T^. „ = we get 



p + d(p + p r ) + n/3(p T + p) = 
p' r + n/3' (p r - p T ) = 
and a Xz (p r - p T ) = -£-p T , i = l,2,...,n 



(5) 



Now for simplicity, we assume p T = and p r = (7 — l)p, < 7 < 2. As a result, conservation equations 
simplify to 



(*) P = Poo(r, X!,x 2 , x n )e-^ a+n ^ 
(it) a Xi = i.e. a is independent of Xi 's 
(Hi) p r = poo{t, xi, x n )e~ nfi 



(6) 



Assuming for simplicity, 



p T = and p r = (7 - l)p, < 7 < 2, (7) 
and integrating the conservation equations we get (choosing the arbitrary integration functions appropriately) 

p(r,t) = Po^- 



and a = a(r, t) (i.e. a is independent of Xi 's) 
From equations (3), (4) and (7) we have the evolution equation of R as 

R 2 = -R 1 - n [poR 1 -~ < C{r)-< + g(r)\ 
n 

where g(r) is an arbitrary integration function and assuming C(r) = ^ + v' to be independent of 't'. 
Thus choosing g[r) = 0, R has the explicit solution 



(8) 



(9) 



t-U= I I , V 7 . [r— - i?— ] 



The hypersurface t = i s (r) describing shell focusing singularity is characterized by R(t s (r), r) = i.e. 



(10) 



t s {r) -U = 



>2n \ C(r)i ^ 
-.r 2 



n + i vpW) 



(11) 
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Now for regularity near the central singularity (i.e., r = 0) we assume 

C(r)=ET=oC j r^ 
Mr)=Er=oftr" +7+j 



(12) 



Then from (11) the time for the central shell focusing singularity is given by 



to 



2n 



Urn t s (r)=t l + T ^^ 



(13) 



The apparent horizon is characterized by R 2 = 1 and if t a h is the time of formation of apparent horizon then 

/2n 



t a h(r) -t = 



^—{B r + 0(r 2 )) - A r q (l + 0{r)) 
7>o 



where, B 



7 Ci 



Co 



EL 

I'o 



(n + 7) 



(14) 



We note that if B < then t a h{r) < t i.e., trapped surface forms earlier than the formation of singularity 
so we can not visualize it (not naked) but if B > then singularity appears before the formation of apparent 
horizon so it is possible to visualize the singularity (naked). It should be mentioned that t a h > to is the 
necessary condition for the visibility of the singularity (locally) while t a h < to is the sufficient condition for 
the formation of a black hole (globally). 

Further, if C\ = = p\ (i.e., a = 0) and q < 2 (q > 1) then from (14) the dominant term is A r q which 
appears in negative sign in the expression for (t h(r) — to)- Hence the collapse always leads to the formation of 
a black hole. Now the restriction q < 2 leads to n + 7 > 4 which implies D{= n + 2) > 6 for 7 = 0, D > 5 for 
7=1 and D > 4 for 7 = 2. Therefore, the validity of cosmic censorship conjecture in 'D' or higher dimension 
depends on the equation of state. 

Secondly, if the arbitrary function g(r) is non-zero then the evolution equation (9) have complicated solution 
of the form 



t-U 



2n 



1 



1 



/ \ , -, 1 r 2 2-Fi[a, -,a+ 1,- 
g(r) n + 1 I 2 



C(r) 



-7r 1 -7, 



1 2F 1 [o,-,a + l, ] 



< 7 < 1 with a = 



and for 1 < 7 < 2 

t-U 



2n C(r)-y/ 2 



po(r) n + 7 



r— 2 F 1 [b,-,b+l,- 



Po(r) 



-R— 2Fl [b,-,b + l,- 



C{r)~<R'i- 1 g(r). 
Po(r) 



with b = 



2(7-1)- 



Thus proceeding as above the time difference between the form of trapped surface and the central singularity 
is given by 



tah — to 



<i K o 1 a + 1, — ] 91 + (3 + n- 2 7 ) ff0 



5o 



2 F 1 [a+l,f,a + 2,-^> 
^ 5o 



+ 0(r 2 



, (2ff /n) 2 - + - 2 



P r 1 , 1 2 .9o - 1 PoC 

, z F x a, -,a+l, ^~ 

(" + l)v/ffo 2 n ,g 



2(1-7) 1 n + l 
• n-1 "Tn — 1 
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for < 7 < 1 
and 

tah - to 



7 /2 



3/2 



^[b, 1,6+1, + M ^[6, i, 6 + 1, - 

2 po Oo 2 po 



7-1 



(" + 7)C 7 
(n + 37 - 2)p 



2 Fi[6 + lj, 6 + 2, -^£]( 7 C lffo po + Cogipo - goC pi) 
2 Po 



n + 7 \ n 



Po 



2 po \ n J 



for 1 < 7 < 2. 



From the above time difference, it is not easy to speculate the end state of collapse (naked singularity or 
black hole) as it involves so many arbitrary (or initial) constants. However, if it so happens that the coefficient 
of r vanishes then the leading order term in r may depend on 7 if the following inequations are satisfied 

2(1-7) n + 1 

1 < -i + < 2, for < 7 < 1 

n — 1 n — 1 



K^±2+ 2 -^l<2, for 1< 7 < 2 
n — 1 n — 1 

In that case as before we always have black hole as the final state of collapse. The above inequalities show 
that the dimension parameter n depends strongly on 7. In particular, for 7 = 1 we have n > 3 i.e., we always 
have black hole solutions for six and higher dimension which is in agreement with result in dust collapse [11]. 
For different equation of state the restriction on n is as follows: 



TABLE-I 



7 





1 


4/3 


2 


n 


6 


4 


5 


7 



The table shows the least value of n for different value of the equation of state parameter 7 in formation of 
black hole. 



III. GEODESICS AND THE NATURE OF SINGULARITY 



As before, using for simplicity g(r) = = f(r) the scale factor R(t, r) has the simple explicit solution (choosing 
the initial time U = 0) 

2 

(15) 



R = 



2^ C(r)T/2 



Following the geodesic analysis of Joshi and Dwivedi [21] for TBL model, let us introduce the functions: 



5 



X=4: 

t _ 2rB' 
5 ~~ B 

v = rQL 

C = 2S 

S ~ r m{n+y-2) 

Q = er v 



r\ _ n + -, 

V - (m -i) ( „ +T) ' 
r 1 



withm> landB(r) = ^p. 

In this section the visibility or non-visibility of the singularity is characterized by examining the possibility 
to have any outgoing null geodesies which are terminated in the past at the central singularity r — 0. Suppose 
this occurs at singularity t = to at which i?(£ ,0) = 0. For convenience , we consider the radial null geodesic, 
given by 

^=e<*=R' + Rv' (17) 
dr 

Now choosing, U = r m , the above geodesic equation can be written as 

W = U{X,u) (18) 



where U(X,u) = (1 - y/ + I V with H = + xln£=*m 



If we approach the singularity at R — 0, u — along the radial null geodesic then the limiting behaviour of 
the function X is given by 



X = Urn X = lim f = lim ^ = Urn U{X,u) =U(X ,0) 

R^Ou^O R->0u->0 i?^0u^0 R^Ou^O 

This is a polynomial equation in X a with explicit form 

X^-\m - -4-] - Q xt + "- 2)/2 + e VCo + v/Co^ n+7)/2 
n + 7 

Here the suffix '0' stands for the value of the variable at r = 0. 

If it has at least one positive real root then it is possible to have a radial null geodesic outgoing from the 
central singularity and the singularity will be naked. The above polynomial equation shows that whenever 'n' 
occurs it is always with 7 in the form n + 7. So if both n and 7 changes keeping n + 7 to be invariant then the 
final fate of the singularity remains same i.e., if we have space-time dimension D = n + 2 and equation of state 
p = (7 — l)p and we have another space-time dimension D = n + 1 and equation of state p = jp then the end 
state of collapse will be same for both space time. 



Co 



n + 7 



Vo 



(20) 



G 



TABLE II 



7 


m 




Co 


Co 


e 






Positive roots (X ) 




















4D 


5D 


6D 


7D 


8D 


10D 


14D 

















.01 


.1 


.21 


.31 


.45 


.62 


4/3 


1 


-6 


.05 


.01 


-5 


.03 


.14 


.25 


.39 


.41 


.52 


.66 


2 












.1 


.21 


.31 


.39 


.45 


.55 


.67 




4/3 
2 


1 


0.1 


1 


1 


1 



13.6 

- 


- 

- 


— 


— 


- 

- 


- 

- 


- 
- 















1.19 


.11 
.96 


.36 
.9 


.54 
.88 


.66 
.87 


.81 
.86 


- 


4/3 


4 





5 


0.1 


4 


.2 
.93 


.43 
.89 


.58 
.88 


.69 
.87 


.76 
.87 


- 


- 


2 












.36 
.9 


.54 
.88 


.66 
.87 


.74 
.87 


.81 
.86 



















1.04 


.9 


.87 


.86 


.86 






4/3 


1 





1 


0.1 


1 


.2 

.89 


.43 
.87 


.59 
.86 


.69 
.86 


.77 

.85 






2 












.36 
.87 


.54 
.86 


.66 
.86 


.74 
.85 









Due to the complicated nature of the above polynomial an exact analytic solution is not possible for Xq. So 
we study the roots by numerical methods. The above table (Table II) shows the dependence of the nature of 
the roots on the variation of the parameters involved. It is clear from the definition of the functions in eq. (16) 
that the parameter m, £o, Co ar e always positive while r]o, Oo may take both positive and negative values. 
However, it is not possible to find any definite conclusion due to the variation of the parameter 6q. 



A. Strength of the naked singularity 

The strength of a singularity (Tipler[22]) is characterized by examining the state of a body that falls in it. 
If all objects that fall into a singularity are destroyed by crushing or tidally stretching to zero volume then 
the singularity is called gravitationally strong otherwise it is called weak singularity. A precise mathematical 
definition was given by Clarke and Krolak [23] . According to them, a sufficient condition for a strong curvature 
singularity is that, for at least one non space- like geodesic with affine parameter A — > on approach to the 
singularity, we must have 

Urn X 2 R ij K i Ki>0 
A^0 [2L) 

with K l = the tangent vector to the radial null geodesic. For future-directed radial null geodesies that 
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originate from the naked singularity the above limit can be written as (using L'Hospital's rule) 

nCo(H -noX )(£o-^n(n+~/)rio) 



Urn X 2 h',.l\- 1\> 

< > (22) 



where H = H(X , 0),N = N(X , 0). 

The singularity is gravitationally strong in the sense of Tipler if 

1 , L (n + 7)O ] 

Co - 2 n ( n + VVo > max < 0, ^ — ^ 

I X 2 ) 



or 



1 , x L (n + 7)6 
so - 2 n ( n + 7)% < mzn < 0, ^ — 



If the above condition is not satisfied for the values of the parameters then \™ § \ 2 RijK l K 3 < and 
the singularity may or may not be Tipler strong. 

IV. DISCUSSIONS AND CONCLUDING REMARKS 

The paper deals with gravitational collapse in (n + 2) dimensional quasi-spherical Szekeres' space-time for 
perfect fluid with barotropic equation of state p r = (7 — l)p. The role of the parameter '7' has been discussed 
in formation of naked sigularity. It is observed that the end state of collapse (black hole or naked singularity) 
remains invariant as long as the sum 'n + 7' is fixed i.e., a perfect fluid with equation of state p r = p in 'n' 
dimension, dust in (n + 1) dimension and exotic matter with equation of state p r = — p in (n + 2) dimension 
are equivalent from the point of view of final state of collapse. This interesting feature has been shown both 
for local and global nature of the singularity. However, one may note that the above results are valid for the 
specific choice of C(r) in the series for given in the eq.(12). Finally, the strength of the naked singularity has 
been examined using the criterion introduced by Tipler [22] . It is found that naked singularity may be a strong 
curvature singularity or not depending on the choice of the parameters at r = 0. 
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